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Introduction 

A  dropbar  system  has  been  developed  to  produce  stress  pulses  above  1  GPa 
with  durations  of  about  500  ps.  These  pulses  are  of  similar  shape  and 
amplitude  to  those  encountered  in  large  explosive  tests.  This  system  was  to 
have  been  used  to  test  a  stress  level  indicator  based  on  the  electronic  phase 
transition  which  occurs  in  cerium  metal;  however,  due  to  the  extreme 
reactivity  of  Ce, no  gauges  could  be  constructed.  We  observed  reactions  with 
epoxy  over  a  period  of  several  days  which  may  explain  the  erratic  results  of 
a  previous  study. 

We  recommend  that  future  work  on  stress  level  indicators  exclude  Ce 
because  of  its  high  reactivity.  Among  the  metals,  T1  ^In  ^  and  BiSn^  both 
have  phase  changes  below  1  GPa  and  warrant  further  study.  Another  inter¬ 
esting  class  of  metal  phase  transitions  are  melting  transitions  many  of  which 
have  pressure  dependence  such  as  Ga,  Hg,  Hg-Tl  and  Rb-Na.  Semi-conductor 
metal  transitions  should  also  be  investigated. 

2.  Background 

Design  of  hardened  underground  structures  for  defense  relies  on 
accurate  constitutive  models  of  the  rocks  or  soil  surrounding  such 
structures.  The  principle  sources  of  data  for  these  models  are  small  labora¬ 
tory  tests  and  large  field  tests.  However,  confidence  in  the  data  thus 
obtained  has  been  slow  to  develop  because  of  uncertainties  in  the  effects  of 
scale  on  rock  properties,  on  the  one  hand,  and  in  the  calibration  of  stress 
gauges  used  in  field  tests,  on  the  other.  This  research  and  development 
project  was  undertaken  to  provide  a  real-time,  in  situ  calibration  technique 
for  stress  gauges.  The  objective  was  to  produce  an  electrical  signal  in  the 
form  of  a  box  car  fuction  whose  upward  and  downward  edges  would  correspond  to 
known  stress  levels.  A  subsidiary  objective  was  to  develop  a  gauge  cali¬ 
bration  facility  where  gauges  to  be  used  in  the  field  could  be  subjected  to 
transient  stress  pulses  of  an  amplitude  and  duration  typical  of  those  en¬ 
countered  in  field  testing. 


The  difficulties  of  credibly  calibrating  piezoresisti ve  stress  gauges 
are  familiar  to  anyone  who  has  used  them.  The  two  most  commonly  used  metals 
are  manganin  and  ytterbium  (Yb).  The  piezoresistive  coefficient  of  both 
metals  depends  on  whether  the  material  has  been  hot -worked  such  as  drawn  wire 
and  hot-rolled  foil,  cold-worked  as  is  the  case  of  most  foils,  or  annealed. 
Furthermore  both  materials  may  exhibit  hysteresis,  a  different  response  to  a 
given  stress  depending  on  whether  the  stress  is  increasing  or  decreasing. 

The  measured  resistance  can  also  be  changed  by  simple  geometric  changes  at 
zero  stress  such  as  stretching  or  puncturing.  Therefore,  a  need  exists  for  a 
technique  to  provide  real  time,  i_n  situ  caliDration  for  such  gauges  when  they 
are  to  be  installed  in  field  tests.  One  would  like,  for  instance,  an 
indicator  based  on  a  phenomenon  that  depends  on  stress  alone  and  not  on 
strains  or  metallurgical  history. 

The  approach  to  the  main  objective  was  one  used  by  Physics  Inter¬ 
national  (PI)  from  1974  to  1977*-1-*.  The  basis  of  their  work  was  the 
existence  of  phase  changes  in  a  wide  variety  of  materials  including  both 
metals  and  insulators.  Literally  hundreds  of  phase  changes  have  been  found 

to  occur  in  a  wide  variety  of  materials  when  they  are  subjected  to  high 
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pressures.  The  work  of  Bridgman*-  J  is  by  far  the  most  extensive  by  a 
single  investigator  in  this  area,  but  numerous  other  workers  can  also  be 

The  properties  that  an  ideal  stress  level  indicator  phase  change  should 
have  can  be  described  as  follows:  (1)  The  phase  change  should  be  rapid, 
reproducible  and  reversible  (but  not  necessarily  at  the  same  pressure), 
resulting  in  a  large,  discontinuous  change  in  a  measureable  physical 
property.  (2)  The  pressure  of  the  phase  change  and  the  magnitude  of  the 
discontinuity  should  be  independent  of  temperature.  (3)  The  material  should 
be  readily  available,  non-hazardous,  and  relatively  stable  under  normal 
conditions  of  atmosphere,  temperature  and  humidity.  Probably  no  material 
meets  all  of  these  requirements  for  ideality;  in  order  to  be  successful  the 
first  criterion  must  be  met,  at  least.  As  we  shall  see,  none  of  the  metals 
investigated  by  PI  fully  meet  the  second  and  third  criteria. 


The  study  of  PI  concentrated  on  two  particular  types  of  phase 
changes  -  transitions  in  metals  exhibiting  a  discontinuity  in  electrical 
resistivity  and  transitions  in  dielectrics  exhibiting  a  discontinuity  in 
specific  volume  and  dielectric  constant  which  could  lead  to  a  capacitance 
change. 

After  surveying  26  materials,  they  selected  four  metals  (Bi,  Ce, 

T1  77In  23.  Bi^Pb),  one  semiconductor  (HgSe),  and  five  insulators  (KNO^, 

NH^I,  ortho-nitrochlorobenzene,  camphor  and  thiourea).  After  a  series  of  gas 
gun  tests,  PI  recommended  Bi,  Bi^Pb,  Ce,  thiourea  and  NH4I  for  further  study. 
These  materials  experience  phase  transitions  at  about  2.6  GPa,  2.5  GPa 
(loading)/!). 5  GPa  (unloading),  700  MPa,  350  MPa  and  50  MPa,  respectively. 

HgSe  and  Tl/In  are  both  toxic  and  the  conductivity  of  the  Tl/In  alloy  is  too 
high  to  make  a  good  gauge. 

Selection  of  a  material  for  further  investigation  from  the  PI  subset 
was  constrained  by  two  further  considerations.  First,  thiourea  and  NH^I  were 
rejected  because  it  was  felt  that  resistivity  measurements  presented  fewer 
experimental  problems  than  capacitance  measurements.  The  subsidiary 
objective  of  developing  a  facility  for  calibration  of  field  gauges  led  to  tne 
rejection  of  the  two  Bi  materials  because  the  high  pressures  required  could 
not  be  obtained  in  a  drop  bar  apparatus.  Given  these  constraints,  cerium  was 
selected  as  the  material  for  further  investigation. 

Cerium  exists  in  a  face  centered  cubic  (a)  phase  at  300K  and  0.1  MPa. 
However,  when  cooled  below  about  263K  it  begins  to  invert  to  a  hexagonal 
close  packed  (g)  structure,  and  below  about  100K  inverts  to  a  second  face 
centered  (a1)  structure  13.5  per  cent  denser  than  the  a  structure.  The 
difference  between  the  a  and  a'  structures  is  attributed  to  an  electronic 
transition  between  the  4f  and  5d  subshells  of  the  Ce  atom  and  is  accompanied 
by  a  marked  decrease  in  resistivity  of  Ce.  None  of  these  transitions  is 
isothermal;  the  ratio  of  a  to  3  or  B  to  a'  varies  continuously  over  a  range 
of  temperature  of  up  to  40K.  Furthermore,  there  is  a  definite  hysteresis 
depending  on  whether  temperature  is  increasing  or  decreasing.  It  is  also 
seen  that  working  above  263K  facilitates  the  6  to  a  transition  and  working 
below  77K  facilitates  both  the  g  to  a'  and  a  to  a'  transitions. 
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Bridgman1"  ’  J  observed  a  phase  transition  and  a  resistivity  change  at 
0.7  GPa  in  Ce,  and  Lawson  and  Tang^  determined  that  the  high  pressure  phase, 
at  1.5  GPa  and  room  temperature,  had  a  face  centered  structure  with  nearly 
the  same  lattice  constant  as  the  a1  phase.  Herman  and  Swenson^  later  made 
volume  measurements  at  various  pressures  and  showed  that  the  high  pressure 
phase  is  continuous  with  the  low  temperature  a'  phase.  A  hysteresis  of 
almost  0.2  GPa  was  observed  between  loading  and  unloading.  The  P-T  phase 
diagram  of  Ce  based  on  their  work  is  shown  in  Figure  l.  The  high  pressure 
extension  of  the  a  -  0  transitions  is  purely  conjectural.  To  further 
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complicate  matters,  Swenson1  J  has  suggested  that,  at  about  1.8  GPa  and  600K, 
there  is  a  critical  point  at  which  the  a  and  a'  phases  are  indistinguishable. 

Based  on  the  static  pressure  measurements,  Ce  hardly  appears  to  be  a 
likely  candidate  for  a  stress  indicator,  especially  in  view  of  the  sluggish, 
continuous  nature  of  the  phase  transitions.  However,  the  fact  that  working 
the  material  seems  to  facilitate  equilibrium  is  perhaps  encouraging.  In 
fact,  the  PI  studies  may  be  interpreted  as  showing  that  the  phase  change  can 
procede  at  least  to  the  extent  of  20  to  30  per  cent  a  converted  to  a'  in  a 
loading  time  of  less  than  0.5  ys.  However,  the  phase  transition  was  not 
observed  in  all  tests  in  which  it  should  have  occurred.  Their  published 
results  are  not  adequate  to  determine  closely  the  pressure  of  the  observed 
transition . 

3.  Experimental  Techniques 
3.1  Dropbar  Facility 

Dropbar  techniques  have  been  used  for  several  years  to  test  and 
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calibrate  a  variety  of  ground  shock  instruments  .  Dropbars  have  been 
reliable  tools  for  producing  controlled  reproducible  stress  profiles  with 
durations  of  0.5  to  12  ms,  durations  similar  to  those  encountered  in  many 
field  applications.  The  shape  of  the  stress  pulse  can  be  predetermined  by 
the  variation  along  the  bars  axis  of  its  cross-section.  The  duration  of  the 
pulse  (T)  is  the  time  required  for  the  stress  pulse  to  propagate  up  and  down 
the  bar  of  length  L, 


where  V  is  the  wave  speed.  In  our  case  the  bars  was  always  below  its  yield 
point  so  V  =  E/p,  the  rod  wave  velocity,  where  E  is  Young’s  modulus  and  p  is 
the  density.  The  shape  of  the  stress  pulse,  o(t),  is  related  to  the  shape  of 
the  bar,  in  terms  of  its  cross-section,  A(x),  by  the  approximate  relation 

o(t)  *  k  C1A(-|  TV)  (1) 

where  C.  is  the  impact  speed  of  the  bar  and  k  =pV/A*  with  A*  being  the  area 
in  the  gauge  plane  which  supports  the  bar.  In  actual  practice  there  is 
dispersion  in  the  wave  propagation  so  that  c(t)  will  be  a  smoothed  version  of 
A(x) . 

In  our  tests  we  planned  to  use  Ce  gauge  elements  about  2  cm  square  so 
that  a  large  A*  would  be  required.  This  necessitated  the  construction  of  a 
larger  dropbar  facility  than  used  in  previous  tests.  Furthermore,  it  was 
desired  to  work  with  peak  stresses  up  to  1.5  GPa  which  necessitated  use  of  a 
hardened  steel  bar  with  a  yield  strength  in  excess  of  that  value  (220  ksi) 
and  an  impact  velocity  of  about  10  ms-1.  One  other  requirement  was  that  tne 
stress  should  pass  through  the  Ce  phase  transition  gradually  (about  lO^GPa 
s’1  or  less)  so  that  the  effect  of  loading  rate  on  the  sharpness  of  the  phase 
transition  could  be  observed. 

Based  on  these  requirements  we  designed  a  dropbar  as  shown  in  Figure  2. 
The  bar  is  a  0.1016  m  diameter  round  bar  of  300  M  (var)  steel,  machined  as 
shown,  with  an  impact  diameter  of  47.9  im,  and  hardened  to  Rockwell  C52 
(yield-  220  ksi).  The  flat  and  flag  on  one  side  were  machined  to  permit  use 
of  a  velocity  detector  described  below. 

In  order  to  achieve  an  impact  velocity  of  10  m/s  using  free  fall  at  10 
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ms  a  height  of  about  5  m  was  required.  The  dropbar  system  was  constructed 
by  modifying  an  existing  30  ft  (9  m)  drop  tower  at  a  nearby  plant.  The  bar 


Figure  2.  Drop  bar  used  to  produce  calibration  pulses 
Dimensions  are  in  millimeters. 
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is  suspended  from  a  cable  by  a  remotely  operated  catch.  The  bar  is  guided  to 
the  impact  point  by  three  pairs  of  trolleys  which  run  on  three  vertical  rods 
which  are  strapped  together  at  frequent  intervals  for  rigidity.  The  rods  are 
attached  to  the  supporting  beam  at  the  top  and  bottom. 

The  velocity  of  the  bar  is  determined  by  measuring  the  time  required 
for  a  flag  of  known  length  (46.8  mm)  to  traverse  a  light  beam.  This  is 
monitored  by  coupled  GaAs  infra-red  emitting  diode  and  silicon  photo¬ 
transistor  (TIL138)  mounted  on  the  lower  support  bar  of  the  guide  rods.  The 
flag  is  1  nm  thick  while  the  gap  of  the  sensor  is  3.1  mm,  permitting  ail  mm 
alignment  error  for  the  bar.  The  velocity  of  the  bar  is  measured  both 
downward  and  upward  on  the  first  rebound.  This  permits  calibration  of  the 
gauge  as  described  in  Section  4.1. 

3.2  Ce  Gauges 

The  Ce  gauges  were  designed  to  be  a  grid  with  a  length  of  90  run  in  an 
2 

area  of  100  rim  .  Our  first  plan  was  to  etch  the  grid  out  of  thin  foils 
mounted  on  a  mylar  substrate  using  procedures  described  by  PI^.  The 
thinnest  commercially  available  foils  were  25 u m  thick  which  would  yield  a 
resistance  of  only  3  ft  if  the  grid  lines  were  0.1  mm  wide.  Although  that 
resistance  was  only  about  one  tenth  that  desired  it  would  have  been  adequate 
for  a  laboratory  test  environment.  However,  the  foils  were  completely 
oxidized  on  receipt.  Foils  were  reordered  at  a  thickness  of  0.1  mm  which 
would  have  yielded  a  value  of  about  0.75ft  ,  a  marginally  useful  gauge.  This 
approach  also  proved  unsuccessful  because  the  foil  reacted  chemically  with 
the  epoxy  used  to  mount  it  on  the  substrate.  The  reaction  rate  was  slow, 
with  about  6  days  being  required  to  completely  alter  the  metal.  This  slow 
reaction  with  epoxy  may  explain  the  variable  results  of  PI  with  Ce  gauges. 
Those  gauges  used  within  a  day  or  so  after  first  gluing  should  give  positive 
results  (shot  15^)  whereas  those  used  after  several  days  storage  should 
give  negative  results  (shots  5  and  6^). 


Because  of  the  difficulties  experienced  with  Ce  foils,  a  second 
production  scheme  was  devised.  We  attempted  to  deposit  Ce  vapor  onto  a  mylar 
substrate  in  a  vacuum.  Numerous  attempts  were  made  with  various  deposition 
rates  but  none  were  successful.  At  high  disposition  rates  the  substrate 
either  became  hot  enough  to  flow  or  reacted  directly  with  the  Ce  upon 
deposition.  At  low  deposition  rates  the  Ce  atoms  reacted  with  the  residual 
gases  in  the  sample  chamber  before  deposition.  In  either  case  the  resulting 
grid  was  not  a  conductor.  Attempts  to  cool  the  substrate  during  deposition 
were  also  unsuccessful. 

3.3  Yb  Gauges 

Our  original  plan  was  to  combine  Ce  gauges  and  Yb  gauges  into  a  single 

gauge  package  as  shown  in  Figure  3.  The  Yb  would  provide  a  measure  of  the 

stress  history  for  comparison  with  the  output  of  the  Ce  gauges.  The  YD 
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gauges  used  were  manufactured  by  S  .  They  were  of  0.05  mm  cold  rolled  foil 
backed  with  0.25  mm  fiberglass  and  had  a  nominal  resistance  of  20ft. 

The  package  finally  used  contained  only  a  single  Yb  gauge  because  we 
were  not  successful  in  producing  Ce  gauges.  The  Yb  gauges  were  laminated 
between  two  2.5  mm  sheets  of  4130  steel  hardened  to  Rockwell  C40.  The 
laminated  package  was  held  together  by  bolts  around  the  edges.  The  resulting 
gauge  package  survived  repeated  loads  up  to  1  GPa  on  the  dropbar. 

3.4  Data  Acquisition  and  Analysis 

3 

The  gauges  were  powered  by  a  pulsed  power  supply  manufactured  by  S  , 
which  effectively  produced  a  constant  current  through  a  100  ft  resistor  in 
series  with  the  gauge.  The  effective  circuit  is  shown  in  Figure  4.  The 
voltage  drop  across  the  gauge  was  read  with  a  Nicolet  Explorer  3  digital 
storage  oscilloscope.  After  each  shot  the  memory  of  the  scope  was  read 
onto  tape  for  later  processing. 
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Voltage  records  were  converted  to  relative  resistance  in  accordance 
with  the  circuit  of  Figure  4, 


AR(t)  .  R(t)-Ro  _  RoCVRo]Cv(t)-vo] 


*0  K0  o 

Rq  is  the  initial  gauge  resistance,  R^  is  the  series  resistance,  VQ  is  the 
initial  rise  in  voltage,  and  V(t)  and  R(t)  are  the  measured  voltage  and  the 
resistance,  respectively,  at  time  t. 

Relative  resistance  records  were  then  converted  to  stress  Dy  two 
techniques.  In  the  first  the  equations  of  Ginsberg  et  al.^  were  used: 

a  =  .1082  [l-exp(-20.8AR/R0)]  +  .9168AR/Rq  (3) 


for  loading  and 


o  =  o 


>k  -  R 


pk  Rpk  -  Rf 


for  unloading  where  apk  is  the  peak  stress  from  equation  (3),  Rpk  is  the 
resistance  at  peak  stress  and  Rf  is 

Rf  =  Rq  +  ARh  =  Rq  (1.1179  +  0.0542  In  o  k).  (4b) 

In  the  second  technique  the  change  in  momentum  of  the  bar  is  related 
to  the  impulse  at  the  gauge  plane,  with  the  further  assumption  that  the  re¬ 
lation  between  AR/Rq  and  o  is  bilinear,  that  is,  linear  on  both  loading  and 
unloading  with  a  different  proportionality  constant. 


:,-+cj  =  k,  [ 


R(t)-R„ 


+  k  of 


R(t)-R, 


where  Cr  is  the  rebound  speed  of  the  bar,  kj  and  k2  are  the  proportionality 
constants  for  loading  and  unloading  respectively,  and  tg,  t  and  t^-  are  the 
times  of  start  of  pulse,  peak  stress  and  end  of  pulse,  respectively.  Using 
kl^Rpk'R0^R0=opk=k2^Rpk*Rf ^Rf  one  of  the  constants  can  be  eliminated  to  get 


n(t)  - 


) 


with 


M<C  +C  ) 

r  f  ' 

Ro 

V*  I 

A*K 

•  '  i  'pk 

R(t)-Rf 

M(Cc+Cr) 

R  . -R 
pk  o 

Ro 

A*K 

R  .  -R,  * 

pk  f 

(5) 


t  >  t  . 
-  Pk 


/' T*  -  S3  /  ^ 


dt 


pk 


The  principle  difficulty  with  equation  (5)  is  that  A*  is  not  accurately  known 

it 

for  the  flat  pack  gauge  design.  Data  obtained  below  are  reported  assuming  A 

-3  ? 

is  equal  to  the  impact  surface,  i.802  X  10  m  . 

4.  Results  and  Oiscussion 

4.1  Dropbar  Tests 

Figure  5  shows  A R/RQ  vs.  time  for  three  typical  dropbar  tests  at  about 
0,3  GPa,  0.5  GPa  and  1.0  GPa.  In  general,  they  have  the  desired  charac¬ 
teristics  of  peak  stress,  duration  and  loading  rate.  Shots  6  and  7  do  not 
exhibit  the  flat  top  expected  (Figure  2)  and  observed  in  shot  4.  This  is 
probably  due  to  excessive  interaction  of  the  bar  with  the  guide  rod  caused  by 
loosening  of  the  trolley  tracks  which  produced  failure  to  trigger  on  shot  8. 
However,  the  stress  pulses  are  certainly  of  the  correct  shape  for  calibration 
of  gauges  to  be  used  in  large  explosive  tests. 
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time  (Msec) 


Relative  resistance  ( AR/R )  as  a  function  of  time  for 
three  typical  tests  with  peak  stresses  of  about  0.3, 
0.5  and  1.0  GPa.  Shot  numbers  are  shown  below  peaks 


\ 


Table  1  compares  values  for  peak  stress  derived  from  equations  (1),  i3) 

and  (5).  The  lack  of  agreement  is  striking.  The  values  of  equation  (3)  are 

★ 

to  be  preferred  because  the  value  of  A  is  uncertain.  Furthermore, 

derivation  of  both  equations  (1)  and  (5)  assume  that  the  stress  is  dis- 

★ 

tributed  uniformly  over  the  area  A  ;  however,  using  the  solution  of 

r13l 

Boussinesq  for  the  stress  in  a  half  space  subjected  to  a  circular  loadL  J  we 
have 

o( r)  »  1/2  a  Pavg/^V 


so  that  a  (r=0)=l /2  Pavg  When  yielding  at  the  edqe  of  the  load  is  taken  into 
account  a (r=0)  will  be  somewhat  higher,  perhaps  between  1/2  Par  and  2/3  Par. 
Thus  the  values  in  (1)  and  (3)  should  be  high  by  about  1-1/2  to  2  times,  as 
they  are.  This  difficulty  could  be  eliminated  by  machining  a  ring  out  of  tne 
gauge  to  eliminate  the  coupling  to  the  sides.  The  resulting  gauge  might  be 
simpler  to  calibrate  but  perhaps  less  durable. 

The  hysteresis  shown  in  all  the  Yb  gauges  is  considerably  less  than 

[  12  "J 

that  observed  by  Ginsberg  et  ar  J  given  by  equation  (4b).  There  are  two 
reasons  for  this.  The  very  stiff  packaging  apparently  reduces  plastic 
working  of  the  gauge  and  hence  lowers  the  hysteresis.  Furthermore,  only  shot 
6  had  a  virgin  gauge.  In  shot  4  the  gauge  had  been  loaded  three  times  to 
about  the  same  level  during  the  preceding  hour.  Shot  7  used  the  same  gauge 
as  shot  6  after  a  delay  of  5  to  10  minutes.  It  has  been  shown  that 
hysteresis  in  Yb  is  virtually  eliminated  during  repeated  loadings^1*^. 

4.2  Ce  Gauges 

We  were  not  successful  in  producing  any  Ce  gauge  due  to  the  extreme 
reactivity  of  the  metal.  We  observed  reaction  with  epoxy  over  a  period  of 
several  days  which  may  explain  why,  in  earlier  studies,  some  gauges  worked 
while  others  did  not. 


Table  1. 

Test 

Peak  stresses 

in  Figure  5  using 

different  equations.'' 

Stress  (GPa) 

eq.(l) 

eq. (3) 

eq. (5) 

4 

0.25 

716 

328 

565 

6 

0.40 

1011 

450 

1050 

7 

0.95 

1510 

987 

1460 

■i 
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